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We present a general model which includes the ferromagnetic Kondo lattice and the 
Hubbard model as special cases. The stability of the ferromagnetic state is investigated 
variationally. We discuss the mechanism of ferromagnetism in metallic nickel, emphasizing 
the importance of orbital degeneracy and the effect of the Hund's-rule coupling. 

§1. Introduction 

The study of itinerant ferromagnetism is one of tiie most difficult problems in 
solid-state theory. In the 1960's, this problem was studied for the single-band Hub- 
bard model by Kanamori,llP Gutzwiller& and Hubbard. & In the Hubbard model, it 
is now widely known that the ferromagnetic ground state is realized only in excep- 
tional situations. An example is the Nagaoka ferromagnetic state, □) which is shown 
to be the exact ground state when a single carrier is doped in the half-filled U = oo 
Hubbard model. The other cases were discussed by Tasaki and Mielke,u) who in- 
vestigated the model in which the relevant band is flat or almost flat. Furthermore, 
Penc et al.i-* discussed ferromagnetism in such a way that their formulation includes 
flat-band limit as a special case. Our theory presented here is also aimed at gen- 
eralizing the model for itinerant ferromagnetism so as to include the Nagaoka state 
in the Hubbard model. Our generalization is made by including the effect of the 
Hund's-rule coupling, rather than by considering the specific lattice structure. 

In a pioneering work on the electron correlation effect, KanamoriEl^ discussed 
ferromagnetism of Ni. Although he started with degenerate 3d bands, he finally 
reduced the problem to that of a single band model. This is partly due to over- 
estimation of the interaction energy, which was regarded as being roughly twice as 
large as the band width. With respect to this point, it is now commonly accepted 
that the interaction energy and the band width in Ni are of the same order of mag- 
nitude. i^'S'llj^ Therefore, it is necessary to reconsider the effect of the Hund's-rule 
coupling. 

In ferrom agn etic nickel, satellite structure was observed at 6 eV below the Fermi 
energy. This structure is interpreted as a two-hole bound state. Illt'llB' To 
reproduce the structure, one must take into account electron-electron and electron- 
hole multiple scattfiring effects. J'or this purpose, formal treatments have been given 
by many authors, ll3)'E3)'llM)'E2P'EI) and the resulting mathematical expressions are all 
in similar forms. In confrmTiity with these results, several works have been presented 
to explain the structure. £3 'E^^'ES The conclusion is that it can be explairied iii the 
single-band model if one assumes the Hubbard repulsion U^W ~4 eV. t3*'ll3)'EfP 
However, it then becomes difficult to maintain the complete ferromagnetic ground 
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state for such a small value of Therefore the effect of degenerate d-orbitals 

and the Hund's-rule coupling should be considered. 

We investigated itinerant ferromagnetism in 3d transition metal systems in the 
generalized Gutzwiller approximation, Ej^ and concluded that the atomic correlation 
effect is quite effective in realizing ferromagnetism in Ni. In this paper, more precise 
discussion on the local stability condition of the complete ferromagnetic state is given. 
We restrict our argument to absolute zero temperature. Not only the longitudinal 
but also transverse component of the Hund's-rule coupling is taken into account. We 
reproduce the theoretical expressions derived for the Hubbard model by considering 
the multiple magnon scattering effect S^'S^'El'' Ell* as a special case of our general 
result. 

In the next section, we present a general model to treat itinerant ferromagnetism. 
The model includes the ferromagnetic Kondo lattice model, which was treated in our 
previous paper. IHS In §§^ and|^, the formulae required in the following sections are 
derived. As an example, we discuss a two-band model, which consists of a pair of 
Hubbard models coupled by the Hund's-rule coupling (§^). Then a simple theory 
for metallic nickel is given in It is concluded that the Hund's-rule coupling effect 
is indispensable. Conclusion is given in the last section, §^. 

§2. General model for the itinerant ferromagnet 

As a general model for an itinerant ferromagnet, we begin with the Hamiltonian 

H = T + V, (2-1) 

t^^r^ c ■ 



y = fii^^fii^l + U' ^ nifj_hiu - J ^ Sifj,- Siu, (2-2) 



where 



^E4/.<7'^<T,<T'Ci^a', (2-3) 



fj.fj' 



with the Pauli matrices cr. The subscripts i, j and fj,, v designate lattice sites and 
localized orbitals, respectively. The summation with < [x ^ v > denotes that the 
sum is taken over pairs of orbitals (//, v). 

For simplicity, we assume that each component of degenerate bands is labeled 
by the index ^, or i^^ = tfjdf^p. Then we have 



T= J2 ^kAf.a^^kf.a, (2-4) 
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where 

£fcM = -E*f.e''^'^"''^- (2-5) 
j 

If we further assume a tight-binding band, or t = t^ for the nearest-neighbor 
pair i,j and zero for any other pairs i ^ j, we have 

ek, = -t,Y.^"''-iu^ (2-6) 

where the summation is taken over the nearest neighbor vector 5. The second term 
is a constant shift of energy. In the following sections, we assume tf^ = and use the 
//-dependent Fermi energy ej^. Then, 



J2^kf. = 0. (2-7) 



k 

This assumption does not cause any difficulty since we are prohibiting the inter-band 
transfer of particles by the assumption made above (2-4). 
The interaction part is rewritten as 



where 

or, explicitly, 



1,11 



As the ferromagnetic ground state itself is weh described by band theory, i3) it is 
safely assumed that the electron correlation effect is irrelevant in the ferromagnetic 
state. In fact, this is true in the single-band Hubbard model, where the Pauli prin- 
ciple makes interaction irrelevant. If the complete ferromagnetic ground state is 
uncorrelated, an analytical treatment to estimate the elementary excitation energy 
becomes feasible. Therefore, for this purpose we assume U' = J/4. Then, 

V = f/^nj^ini^i 

J2 + SiM^^v - (^iA^T^M + ni^,lh^u^)) . (2-11) 
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Finally we have two positive parameters U and J to describe the interaction 
effect. In our model, the term energy of the singlet state formed by two electrons in 
the same (different) orbital is U (J), while that of the triplet state is zero.Es^ For 
example, for metallic nickel, the relative order is J'^U'^W ~ 4 eV, where W is the 
band width. In the following sections, we investigate the model Hamiltonian which 
consists of ( ^-41) and (2-11), i.e., 

H = T + V 

- ^ {4f,Siu + Si^4u - (^*mT^M + ^iMi^ii^t)) • (2-12) 

The complete ferromagnetic ground state |F) is given by 

n 4tIO)' (2-13) 

where |0) denotes the vacuum and Ef^ is the Fermi energy. 

The above model contains several models which have been intensively discussed 
in the field of the itinerant ferromagnetism. For example, we recover the single-band 
Hubbard Hamiltonian by neglecting orbital degeneracy. On the other hand, we 
reproduce the ferromagnetic Kondo lattice model with the localized spin Sj = D/2. 
For this, we may assume = 1 for bands with fi = 1, 2, 3, • • • , -D, and = n for 

the conduction band (0 < n < 1). The completely filled band (n^ = 1) represents 
the Mott insulator and acts as an array of localized spins, = 1/2. These localized 
spins are coupled to form Sf = D/2 hy the on-site Hund's rule. To describe metallic 
nickel, one may assume triply-degenerate bands with = 0.2 for /z = 1,2 and 3 

(§!)• 

§3. Spin wave 

In this and the next section we derive formulae required in the following sections. 
These are generalization of our previous results. 

3.1. Approximation I 

First we consider the trial state for the spin-wave-excited state 

%, = A-'/^Y.^''"'S-\F), (3-1) 

i 

where 

•^r^E^r^- (3-2) 

and A is the total number of lattice sites. Description with this trial state corresponds 
to the random phase approximation (RPA) in the strong coupling limit. 
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The energy expectation value of ( p-lD gives the magnon dispersion, i) 

1 1 

^Oq = ■ -^^nkij^{ek+qfj. - ek^), (3-3) 

k,fi 

where 

nk^^ = (F|c^Cfc^||F), 

k 

2S = J2n^. (3-4) 

In these expressions, re^^ is the step function, is a carrier density in the band ^, 
and 5 represents the total spontaneous magnetization. 
In the long wavelength limit, we have 

LOog = Doq^ = DQ^q\ {q ^ 0) (3-5) 

where 

= (3-6) 

and ^ 

^gti = ^'Yf^k^Ji£k^l■ (3-7) 

A: 

Here we have assumed a tight-binding band with the coordination number z. The 
negative quantity e^^ is the total kinetic energy per site of the band //. Note that 
energy of the ground state |F) is given by 

E^ = AY,eg^ <0. (3-8) 
Thus the magnon stiffness constant Dq is proportional to \E^\: 

°» = M^ 

3.2. Approximation II 

It is well known that the RPA overestimates tha stability of a ferromagnetic 
state. To improve upon (|3-1|) , we postulate the statecB) 

i 

where 

\i) = A-^l^ Y.^!;{r, - r,)clc,, + f-{r, - r,)Q^4)|F), (3-11) 
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where we expressed q^i simply as Cj^. The state is constructed so as not to 
suffer any energy loss by U nor J. Therefore it is used as a trial state for the case 
U = J = oo. The functions fj^{rj —Ti) are determined variationally. Note that (|3-l| ) 
is recovered by assuming fj^{rj — Vi) = dr^-n- 

We introduce the magnon creation operator S]^ defined by 

%^S\\F). (3-12) 

Then the energy of the variational state ( p-lOD is given by 

{nS,[H,S\]\F) 
U)„ = r^^ . (3-13) 

{¥\S,S\\¥) 

After some calculations, we obtain 

(F|5,5t|F) = ^Y. ^lAfk, + E ^kp^fk^fp,, (3-14) 

{F\Sq[H, 5]]|F) = ^ E ^kqp.fkpfkp + ^ E ^kpqfj.fkfifpiJ.^ (3-15) 
k,fj, k,p,fj, 

where fk = n^fk + (1 — n'k)fk = fk + is & sum of the Fourier transform of 
fj^{rj —Ti). Quantities appearing in the above expressions are calculated as follows; 

Al^ = {2S - n^)n^(l - n^^) + (25 -n^ + 1)(1 - n^)nfc^, (3-16) 
^kpfji = (25" - n^)(l - nkp){l - Up^) + {2S -n^ + l)nk^np^, (3-17) 

Akq^ = {2S - n^) {^f,£k^ ~'A^ ^fc'M^fc'M j (1 - ^kfi) 

+ {2S - + 1) ^(1 - nk'^)ek'^ - (1 - n^)ek^ 

+ 1 (1 — n^)'^ek+qfj. — -|2 E ''^k' ^nk" fj,£k' -k" -k-qfi j ra/c/i, (3-18) 
\ k',k" / 

rkpq/i = 2(25 - n^)ekp. (ukfii^ - Up^) - (1 - nk^)np^ (3-19) 

+ E "''^''M i^k'+qfi + ^k'-k-p-qfi) + (1 " n^){ep+q^ + Ek+qp)^ rik^Up^. 

If we assume t^- = for the nearest-neighbor pair i, j, Akq^ and /fcpg/x are expressed 
in terms of e^^, (p-7[), as 



+ (25 -n^ + l)(|eg^| - (1 - n^)£k^nkp, 

+ (^(1 - n^f - i^f^ Sk+q^uk^, (3-20) 
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rkpqpi = 2(25' - n^)ekpi{nkti.{'^ - np^) - (1 - nkf,)np^^ 



+ + ^k+p+qfj.) + (1 - ny){ep+q^ + Ek+q^,) j nkfiUp^j,. (3-21) 

These are Generalization of our previous results for the ferromagnetic Kondo 
lattice model. Ej) Their correspondence becomes clear by introducing the 'localized- 
spin' Sf, 

2Sf = 2S-n^. (3-22) 

The quantity Sf represents the total spin corresponding to all bands other than that 
containing the particle under consideration, i.e., all bands except for /i. Generally, 
Sf depends on the band index /i and can take an arbitrary. value, although it must 
be a half integer in the ferromagnetic Kondo lattice model. c3) We use this parameter 
Sf in the following sections. As noted previously, ll' the case Sf = can be used to 
investigate the U = oo Hubbard model. 

We can estimate the spin-wave stiffness D analytically from ( p-13 ).cJ^ The result 

is 

iOg = Dq^ {q ^ 0) 

= ^Pom - ^D^)q^, (3-23) 

where Dq^ is defined in ( |3-6| ), and 



1 (l-n„)^I„ 



2Sl + {\egfj,\/zt^){If,/2tf,) 



1 ^ 



I, = jE^^k„ (3-25) 

k 



kOiM 



-k, = (3-26) 
These expressions for D, below (3-23), will be used in the following sections. 



§4. Individual-particle excitation 



To describe the individual-particle excitation which has a spin component an- 
tiparallel to the spontaneous magnetization, we introduce the following creation 
operator for the quasiparticle in the band /i: 

cli = ^ E (hvcU + E f^-^,^^4^ ■ (4-1) 

The expression ( [4-l| ) for /lu = and = 5ij appears as a prefactor of c^^ in iP'g. 

In this case, the trial state does not suffer energy loss due to U nor J. To include 
the effect of finite interactions, we introduce a variational parameter h\j and assume 
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the dependence of h'j_^^ on fj, — u. The dependence on rj — ri is further assumed 
to afford a better variational description. This corresponds to ahowing an up-spin 
particle (c^^ in (^^)) to be created at a site different than the site i where a down- 
spin particle is created by s~ . As will be seen below, the effect of multiple magnon 
scattering is taken into account by this generalization. 
The energy expectation value of (|4-lD is given by 



for which we obtain 



(4-2) 



(4-3) 



(F|C,^i[i^,<J|F) 



A 



efc2 



k2 



+ 



Efc/t + S'Atl^ul 



+ 2 ^ 



n. 



>'kx^l 



f^kifi] 



(4-4) 



where /i^^ is the Fourier transform of /ij'^ multiplied by 1 — n^,^. The parameter 
is defined by 

Qfi = Urif, + JSf, 

and represents the Hartree-Fock exchange splitting of our model. 

We remark that we must assume (i) /lu = for [/ = oo, and (ii) /i^^^ = h^^^ 
and h\j = for J = oo, as noted above. This is because the terms depending on 
and/or J in (4-4), which are positive definite, should vanish in the strong coupling 
limit. 



(4-5) 



4.1. Case J ^ oo 

As the above expressions are too general and complicated, we show simpler 
expressions which are used in the following investigation. To investigate the effect 
of a finite Hund's-rule coupling J, we introduce a new variational parameter hj by 



(4-6) 



The parameter hj measures the effect of the Hund's-rule coupling. In particular, 
/ij = for J = 0, and hj = 1 for J = oo. 

All terms but the last term of (4-4) are independent of the phase of the complex 
quantity /ij. Therefore, it is always energetically unfavorable to have /i^^^ out of 
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phase with /i^^^- Thus, we may assume hj as a real quantity, < /ijjj^^ < 1. This 
imphes that the spins Sj^ at site i should couple in phase with each other by the 
Hund's-rule coupling. For /ij = 1, the spin lowering operator ^u^'j-i^i^^u (|4-l|) 
reduces to a quantity proportional to (3-2). 

After taking the functional derivative of ( ^-21) with respect to /lu and /i^^^ and 
eliminating these parameters, one obtains the eigenequation for the excitation energy 



Z:^{k,Ek^l), 
9ti {^-j^lp) 
1 - {9ti+U)P 



(4-7) 
(4-8) 



where 



1 - ?^fcu 



^ k^ (25'//iJ + nf,){Ekf,i - uj'f,k-k^ - ^ku 



n 



1 1 



and 



3^ ^ JSf{l - h,f. 
Noting the identity 

-T X! ^k2v{hi + (1 - /ij) = '^Sfhi + n^, 



(4-9) 
(4-10) 

(4.11) 
(4-12) 



k2,i' 



and comparing uj' with (p-3D, we observe that uj' is the generalization of the 



A"? 

magnon dispersion in the presence of the inter-band spin-spin coupling. For the 
complete coupling /ij = 1, (|3-3|) is recovered, while /ij = gives the single-band 
result for the band /i. 

Since further variation with respect to /ij leads to a complicated expression, we 
leave /ij as a parameter. Thus, the equation (4-7) is to be solved for E^^i as a 
function of /ij, which in turn should be determined so as to minimize Ej.^\^. The 
expressions following (4-7) are conclusions of this subsection. We note that the 
derivation of (4-7) from ( [4TD is exact. While an assumption is made in ( |4-6| ), this 
does not conflict with the variational principle. Introduction of the single parameter 
/ij to describe the Hund's-rule coupling is largely for the sake of simplicity. 

4.2. Case hj = 1 

The case /ij = 1 is realized for J = oo, where the Hubbard repulsion U is always 
ineffective. This is because the strong Hund's-rule coupling forbids two particles to 
stay on the same site. However, there is a situation in which the assumption /ij = 1 
and /lu 7^ is adequate even for a finite interaction J <SQ- We observed that 
this is in fact the case for the double exchange ferromagnet.cB* By this assumption, 
the number of variational parameters decreases by one, and the calculation becomes 
relatively easy. 
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In this case, from the results of the previous subsection, we obtain an eigenequa- 
tion of the form 

Ek^l - Ekf, - S^{k, Ek^i) = 0, (4-13) 

where 

^,ik,u^)= . J' , (4-14) 

and 

^',{k,u^) = ■ (4-15) 



The magnon energy uiQq is defined in 

In these forms, it is possible to regard Z'^(/c,a;) as the self-energy given as a 
result of the multiple scattering of quasiparticle e^^^ off the spin wave cuQk-ki- In 
fact, a particular case of this result, for the single-band Hii^blmrdmi^ 2S = n and 
g = Un, has already been obtained by many authors. lnP'Ej)'E3)'EJ'E2P Although the 
formulae we derived have the same forms as those derived in the other methods, 
it contains a meaning more significant than that of mere re-derivation. Since our 
results are based on the variational principle, we can draw a definite conclusion on 
the instability of the ferromagnetic state. 

§5. Two-band model 

In this section we consider a simple case to show general features of our model. 
We investigate a special case of (2-12), i.e., a model comprising two bands in a 
tight-binding square lattice structure. Then the band dispersion is given by e^^ = 
—2t{cos{kx) + cos(fcj;)) = 1, 2, f > 0), and the band-width is 2zt with z = 4. This 
model includes the Hubbard model and the Sj = 1/2 ferromagnetic Kondo lattice 
model as special cases. Therefore we can treat them in a unified way. Parameters 
used in the following discussion are the interaction U, J, and the carrier density 
of the band /x. The purpose of this section is to investigate the effect of the orbital 
degeneracy and the Hund's-rule coupling. 

5.1. Spin wave 

In this and the next subsection, we consider the strong coupling limit U = J= 00, 
where the formulae given in §|3| are of direct use. In Fig. |, we show the spin-wave 
stiffness constant D{ni,n2) as well as DQ{ni,n2) as a function of ni and n2, where 
Dq is given by ( |3-S| ) and D is defined below (3-23). In Fig. 2, they are shown along 
the symmetry axes in the ni-n2 plane. 

If the effect of band degeneracy is simply additive, one should expect D{n,n) 
= D{n,0) = D{0,n). In fact, this holds for the result Dq in the RPA, for which 
the curve determined by Do(ni,n2)=const becomes convex in the direction (ni,n2) 
= (1, 1) (Fig. 1^). However, this is not the case for the improved estimate D{ni,n2); 
the curve -D(ni, n2)=0 is concave instead, as shown in the dash-dotted contour line on 
the left of Fig. ||. This is due to a singularity present in D{ni,n2) at (ni, 7^2) = (0, 0). 
Here we note that D{n\,n2) = defines the threshold for ferromagnetism. In our 
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Fig. 2. Spin- wave stiffness constant D and Do 
as a function of (ni,n2). 



Fig. 3. Phase l)oundary determined by the 
condition Ekii ~ £f and Ek2i = £/ with 
k = 0. In the region enclosed by the curves, 
the ferromagnetic state is absolutely unsta- 
ble. 



treatment, the result D < represents the absolute instability of the ferromagnetic 
state. Prom the displayed results, we understand that a multi-band model favors 
ferromagnetism more than a single-band model and that one should go beyond the 
RPA to show this multi-band effect. In the RPA, we cannot prove instability at 
all, since Dq > for any (ni,n2) (see Fig. 2). In particular, one should be careful 
in treating the low-density limit owing to the presence of the singularity mentioned 
above. This point is discussed later. 

Next, we note that the line along (n, 0) or (0, n) (0 < n < 1) corresponds to 
the U=oo Hubbard model and the line along (n, 1) or (1, n) to the ferromagnetic 
Kondo lattice model with Sf =1/2 and J=oo. These two cases were discussed in 
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Ref. 26). As a new result, in the direction {n,n), i.e., in the case rii = n2, the 
singular behavior around (0,0) makes the threshold value so small as nc^O.l. (See 
D in Figs. || and 2). We note that the ferromagnetic state is unstable for n < Uc, 
while we may assume n^i ^ 0.2 for Ni. Thus one may simply consider that the fact 
nNi > ric is relevant to the itinerant ferromagnetism in Ni. This point is addressed in 
the next subsection, since the more stringent condition is obtained by investigating 
the individual-particle excitation. 

5.2. Individual-particle excitation 

On the stability of the ferromagnetic state, investigation of the individual- 



particle excitation proc 
wave stiffness constant. 



ces a more stringent condition than that given by the spin- 
V To investigate the case U = J = oo, we use the expressions 
given in §4.2. The eigenenergy E^^i, regarded as a real quantity, is calculated as a 
solution of ( [4-13 ) ^ ( [4-15 ). We determined the critical boundary in the ni-n2 plane 



using the condition E^^^ = ej^ (/i = 1, 2), where we set /c = for the bottom of the 
dispersion e^^. The result is shown in Fig. 3. The Fermi energy £f^ is fixed for a 
given n^. The ferromagnetic state is absolutely unstable when -Efe^j < ^Z/^- 

As mentioned, threshold curves thus obtained give a more stringent condition 
than that determined by the relation D{ni,n2) = 0. (Compare Fig. 3 with the 
contour line D = of Fig. Q.) In particular, in the direction {n,n), the critical 
concentration is estimated as Uc ~ 0.12. In this respect, let us remark again that 
the stiffness D has a singularity, as noted in the previous subsection. This is due 
to the denominator 25 — > as n ^ in the expression for the stiffness constant 
D, ( p-24| ). On the other hand, Do(^i)^2) goes to a finite value in this limit (ni,n2) 
— > (0, 0) (Fig. 2), because the numerator of ( |3-9| ) also vanishes in this limit. Thus, the 
correction AD/Dq is largest in the vicinity of the origin, as seen in Figs. || and 2. On 
the other side, the equation (|4-13D for E'fc^i contains the spin- wave energy a;o(^) in 
the denominator of U^{k,uj), ( [4-15| ). Our calculation based on (44) is exact, so that 



our result gives the exact lower bound for Uc- Nonetheless, in a physical (i.e. non- 
rigorous) sense for improving the approximation, one should replace uJo{q) ~ D^q'^ in 
( 4-15| ) with uj{q) — Dq^. In this replacement, a more stringent andrealistic condition 



will be obtained, although it then loses variational significance.® As a result, we 
expect that the critical boundary in the region where D <C -Do is drastically modified 
as the approximation is improved. For example, from Fig. 2 along the line (l,n), 
i.e., for the double exchange ferromagnet with 5*^=1/2, the correction will be less 
prominent than the case along (n, 0), the U=oo Hubbard model. Around the origin, 
D differs considerably from Dq (Figs. || and 2), so that the critical in a proper 
treatment will take a larger value than Uc ~ 0.12 of Fig. 3. However, it is still 
legitimate to assert that the doubly-degenerate-band model is more favorable to 
ferromagnetism than the single-band model, although in the following section we 
stress that ferromagnetism of Ni cannot be explained only by the effect of orbital 
degeneracy, but a peculiarity of the density of states is required in addition (^S). 



In any case, we note that the correction AD/Dq (3-24) is a quantity of order 
~ 0{1/S) for large 5, since A^q^ in (|3-25|) is proportional to S. Therefore, 
the RPA is quite reliable for models with highly degenerate orbitals and/or a large 
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spontaneous magnetization, as is expected intuitively. 

5.3. The effect of the Hund's-rule coupling 

In this subsection, we assume U to be infinite but J to be finite in order to 
consider the effect of the Hund's-rule coupling. As in the above subsections, the 
model consists of the two Hubbard models in the square lattice, coupled to each 
other by the Hund's-rule coupling. 

We are interested in the effect of the interaction J which couples the two bands. 
Previously it was concluded that the ferromagnetic state in the single-band Hubbard 
model (the limit J ^ of the present model) is unstable beyond some critical hole 
concentration (5c = 1 — Uc, while it is stable for an arbitrary concentration in the 
Sf=l/2 Kondo lattice model (the limit J — > oo). On the other hand, as we saw 
in the previous subsections, the model which assumes ni = n2 is more stable than 
the single-band model, although it would be less stable than the double exchange 
ferromagnet with Sf =1/2. To interpolate these situations in terms of J, the critical 
carrier density determined by the quasiparticle energy is given below as a function 
of J. To this end, we use the formulae given in §4.1. 

The eigenenergy E^^i is obtained as a solution of (4-7) with 

Ui,{k,Ekf,i) = jf,-^, (5-1) 
since 5(^ = 00 for U = 00. Here, p and are functions of hj and (or 5/), and are 



defined in (|4-9| ) and (441), respectively. As noted there, the variational parameter 
hj has been left as a parameter to be determined to minimize the energy. However, 
in determining the phase boundary, minimization on hj is not necessary; we only 
have to calculate J to satisfy E^f^i = £/ for a fixed carrier density n and for various 
values of hj. The envelope formed by a set of curves for every hj then gives the 
threshold J as a function of n, since the variational principle implies instability in 
all cases if it is the case for any of hj. 

We give results for the model (i) for which we set ui = 1 and n2 = n (Fig. 4), 
and for the model (ii) with ni = n2 = n (Fig. 5). In the former, the single-band 
Hubbard model is interpolated to the Sf =1/2 ferromagnetic Kondo lattice model. In 
the latter, it is interpolated as a function of J to the case of two equivalent Hubbard 
models strongly coupled via Hund's rule. 

For the model (i) (see Fig. 4) the region in which the ferromagnetic state is stable 
is quite large. In particular, the Hund's-rule coupling becomes effective even for J as 
small as 0.5zt, a quarter of the band-width W = 2zt. For large J the ferromagnetic 
phase becomes stable for arbitrary n, as expected for the double exchange model. 
The critical boundary (envelope, thick curve) in Fig. 4 is determined by curves for 
relatively small values /ij^O.4, as expected for the boundary in the weak-coupling 
region; curves with larger /ij are completely surrounded by the envelope and do not 
contribute to it. 

For the model (ii) (see Fig. 5) the envelope strongly bends upward around 
J ~ 2zt = W. This implies that the Hund's-rule coupling becomes effective when 
the coupling energy becomes comparable with the band width. The critical hole 
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Fig. 4. Thick line: Threshold for the stability 
of the ferromagnetic state. The model con- 
sists of the U = oo Hubbard model and an 
array of localized spins Sf = 1/2 coupled 
to each other by J. 



Fig. 5. Thick line: Threshold for the stabil- 
ity of the ferromagnetic state. The model 
consists of two equivalent U — oo Hubbard 
models coupled by J. 



concentration 5c increases from 6c ~ 0.4 for the U=oo Hubbard model to 6c ~ 0.9 as 
J increases from zero to infinity. This, and especially a steep increase in 6c around 
J ^ clearly indicate the importance of the inter-band spin-spin coupling. Thus 
it is concluded that the stability condition based on a single-band model becomes 
more stringent than what should be obtained for realistic bands with orbital de- 
generacy. The Hund's-rule coupling effect, if not complete (or even if hj ^ 1), can 
change the criterion for the ferromagnetic instability to a degree that it cannot be 
neglected in investigating real materials. In the next section, we discuss the itiner- 
ant ferromagnetism in Ni, where the importance of the Hund's-rule coupling is again 
stressed. 



§6. Ferromagnetism in nickel metal 



Discussion given above has been based on the expressions derived variationally 
in §§^ and ^ Thus the results indicating instability of the ferromagnetic state are 
exact. However, in discussing itinerant ferromagnetism of metallic nickel, we prefer 
to give a physical argument since a model for Ni has to be simple enough to be 
tractable. A model density-of-states that we adopt below has some characteristic 
features which are regarded to be relevant to ferromagnetism in Ni. We do not 
attempt to estimate the magnon dispersion, for example, by using our expression 
for D, since we are not sure that the model given below is so quantitative as to 
reproduce an appropriate estimate even for the spin wave spectrum. Therefore, we 
regard the spin-wave stiffness constant D as a given quantity and substitute the 
dispersion ujq = Dq^ for in (|4-9| ). The model and parameters used below are set 
so as to mimic metallic nickel. 
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6.1. Simple model simulating Ni 

To describe metallic nickel, we adopt a simple model as a matter of convenience. 
As one of important characteristic, the density of states of Ni at the Fermi level pf 
takes a quite large value. This is doubtless an important factor for ferromagnetism, 
as indicated from the Stoner criterion pjU > 1. Thus, as an example of models 
simulating Ni, we adopt the model due to Edwards and Hertz, EIP i.e., with the 
isotropic dispersion depending on k quadratically: 

r r 

e.={ (6-1) 



k 



27712 



C , ki<k<kr,. 



Here, the constant C is chosen to make Sk continuous at ki. The quantity kj) 
= (Gtt^)^/^ is a Debye-like cutoff chosen so that the band contains one electron state 
per spin per orbital. The parameters mi and m2 are determined so as to meet the 
following conditions, (i) For the Fermi wave- vector kf = ki, the total density per 
band is n = 0.2. (ii) The bandwidth is 4 eV. For the Fermi energy, we introduce 
a parameter ei(= £ki) which characterizes a peculiarity of the density of states. 
The parameter ki in ( |6-1[ ) is uniquely determined by the condition (i). Then, mi 
is set by = ei for a given ei. Parameters ni2 and C are determined to give the 
total band width 4 eV and £ki-jx = £fci+o- If we assume £i = 0.3 eV to mimic Ni, 
then mi = 8.66 and m2 = 1.35.0^ These choices are based on the hole pictureEP in 
which 0.6 d-hole per nickel atom occupies three sub-bands. Moreover we assume the 
quadratic dispersion for the magnon,EJ^ 



= Dq', (6-2) 
where the stiffness constant is set such that Dkf = 0.2 eV. The density of states for 



the model (64) is shown in Fig. 6. 

We investigate a model consisting of three equivalent sub-bands, each of which 
has the dispersion as given above. The Hund's-rule coupling J increases the 
mean- field exchange splitting g defined in ( [4-5D . This effect, however, is due to 
the longitudinal component of the coupling. Therefore, to investigate the transverse 
coupling effect, which has been discussed only little, we take g (or g/n) as a parameter 
for the one-particle bare interaction energy. Thus, for given g/n, J and n, a set of 
equations (4-7) ~ ([4-9D are solved for -Bfc^j^. Results for a single-band model are 
obtained by setting J = and hj = 0. Note that the parameter U = g/n represents 
U = U for the single-band model and U = U+J for the model with triply-degenerate 
orbitals. The parameters used below are g/n, n, J and £i. 

6.2. Results 

First, we show the phase boundary determined by Ek=Q^i = £/ for the degenerate- 
band model, as well as that of the single-band model (the case J = and /ij = 0). 
The critical interaction Uc = gc/n (eV) as a function of the carrier density n is 
shown in Fig. 7 for ei = 0.3 eV. Since we should have g/n'^A eV and n = 0.2 for 
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CO (eV) 



Fig. 6. Bare density of states of our model for 
metallic nickel. In the ferromagnetic state, 
a carrier of total density n = 0.2 per band 
occupies up to the Fermi energy ei (eV). 



Fig. 7. The critical value of g/n as a function 
of n for El = 0.3 eV. In the region below 
the curve, the complete ferromagnetic state 
is unstable. The curve for J = and hj = 
represents the result for the single-band 
model. 




Fig. 8. The critical value of g/n as a, function 
of £1 for n = 0.2. In the region below the 
curve, the complete ferromagnetic state is 
absolutely unstable. The parameters h,j 
shown for J = 1.5 and J = 3 are those 
which give the most stringent criterion. 



Fig. 9. Density of states p{uj) for g/n = 3.5, 
J = 1.5, and hj = 0.2. The carrier den- 
sity is n = 0.2 for each of triply-degenerate 
bands. The dashed curve represents the 
Hartree-Fock result. The Fermi energy is 
£/=£!= 0.3 eV. 



Ni, we conclude that the Hund's-rule couphng is necessary to stabihze the ferro- 
magnetic ground state. Here we give some examples of the results: For n = 0.2 
and g/n = A eV, the ferromagnetic state becomes unstable for J^O.3 eV, and hj 
becomes non-zero for J^l.3 eV, while hj = 0.7 for J = 3 eV. 

After the importance of the the Hund's-rule coupling is shown, the next question 
is that concerning a peculiarity in the bare density of states. We know that a large 
density of states is not necessary for ferromagnetism in the ferromagnetic Kondo 
lattice model. However, in metallic Ni, the carrier density is so low that it is not 



Hund's-Rule Coupling Effect in Itinerant Ferromagnetism 



17 



clear if the Hund's-rule coupling alone can stabilize the ferromagnetic state. To see 
the effect, one may change the parameter ei, which is set above as 0.3 eV. We display 
the phase boundary g/n as a function of £i for n = 0.2 (Fig. 8). Note that states 
of total weight n = 0.2 occupy up to ei (eV) within the band of 4eV wide. Thus a 
small ei represents a strong peculiarity, favoring ferromagnetism. (See Fig. 6.) We 
see in Fig. 8 that the Hund's-rule coupling is necessary to keep the relevant region 
around {ei,g/n) ~ (0.3,4.0) eV ferromagnetic. Still, a large bare density-of-states 
is also required, as concluded in Ref. [25| ). Note that we require ei to be as large as 
4(A;i/A;d)^ = 1.37 eV to give a simple parabolic band, for which mi = m2- Even with 
ei ~ 0.5 eV, the peak structure in the density of states is pronounced, as shown in 
Fig. 6. 

Although the phase boundary in Fig. 7 around n ~ 0.2 does not depend strongly 
on hj, the coupling in terms of hj plays an important role. This point is clearly seen 
in the renormalized density-of-states curve p{io). To investigate whether we obtain 
the satellite structure, we calculate p{uj) defined by 

p{u;) = --Im^Yl — — -, (6-3) 

where E{k,uj) is given by (4-8). In order to estimate the position of a satellite 
structure, one may neglect the momentum dependence of I]{k,u;), or ujg may be 



replaced with the quantity averaged over q: cOg = |Z)/cq.IHb) Then, the calculation is 
straightforward, and the result is shown in Fig. 9 for ei = 0.3 eV. 

It is important to note that we can assume parameters so as to reproduce the 6 
eV satellite structure while preventing the instability of the ferromagnetic state. In 
the figure, we set g/n = 3.5 eV, J = 1.5 eV and n = 0.2. The parameter hj = 0.2 is 
determined to minimize Ek=o^[ ■ These values for g/n and J are physically reasonable 
for metallic nickel. In the calculation with the other set of parameters, we observed 
that the position of the satellite depends on hj, if it is regarded formally as a free 
parameter. For example, for the same set of parameters except for hj, the satellite 
appears at 5 eV and 7.7 eV for hj = and 1, respectively. The weight of the satellite 
also depends on hj and is an increasing function of hj. This is physically interpreted 
as follows: The weight in the low-lying part of spectrum (the part other than that 
due to the satellite) decreases as the local constraint from the coupling J becomes 



effective, or as hj increases. 
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§7. Conclusion 



We investigated the effect of the Hund's-rule coupling and the orbital degener- 
acy on the basis of the variational treatment of a general model for the itinerant 
ferromagnet. Our model includes the ferromagnetic Kondo lattice model and the 
Hubbard model as special cases. It was also applied to the ferromagnetism in Ni. 
To investigate instability due to individual particle excitation, we took into account 
the multiple magnon scattering effect and reproduced some of the results which 
were derived in the Hubbard model by many authors. Generally, in the theory of 
itinerant ferromagnetism, it is important to meet the requirement imposed by the 
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rotational symmetry of the model Hamiltonian. In terms of the diagram technique, 
one should take inki aixoimt the vertex correction in accordance with the Ward- 
Takahashi identity. E^'EHP'O) In this respect, we avoided complexity by making use 
of proper variational states. We derived the required expressions variationally with- 
out any ambiguous assumption. Variational derivation is advantageous because we 
do not have to be concerned with the problem of selecting the type of diagrams to 
be summed over and we are free from a convergence problem of the summed infinite 
series. Moreover, the result for the instability has an exact significance, since it 
follows the variational principle. 

We clarified the physics introduced by the Hund's-rule coupling by treating 
a two-band model and a model simulating metallic nickel. For the former, we 
noted that attention should be paid to the treatment of the low-density limit of 
the degenerate-band model, owing to the singular behavior of the spin-wave stiffness 
constant D. We displayed figures for D which show that the correction to the result 
of the RPA becomes smaller when the spontaneous magnetization becomes larger. 
Investigation of the individual-particle excitation gives a more stringent condition for 
the instability of the ferromagnetic state than that given by the spin- wave instability. 
However, in both cases, we observed a common tendency that a two-band model is 
quite stable compared with a single-band model. As a function of the Hund's-rule 
coupling J, we investigated how two bands are ferromagnetically coupled. In the case 
where one of them constitutes localized spins, J is quite effective, as expected. In the 
other case, where the two bands are equivalent, J stabilizes the ferromagnetic state 
in the regime J^W. We conclude for Ni that a peculiarity of the density-of-state 
curve at the Fermi level is indispensable, while the effect of Hund's-rule coupling is 
also required. In this respect, ferromagnetism in Ni differs from the double exchange 
ferromagnet, in which the complete ferromagnetic state is made stable only by the 
Hund's rule coupling. However, it is conceptually simple to discuss itinerant ferro- 
magnetism on the basis of a unified model as we proposed in this article. It was also 
shown that the effect of the Hund's-rule coupling including its transverse component 
can explain the 6 eV satellite structure as well as the ferromagnetic ground state of 
metallic nickel. 
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